A Note on Tiling under Tomographic Constraints 
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Abstract 



Given a tiling of a 2D grid with several types of tiles, we can count for every 
row and column how many tiles of each type it intersects. These numbers are called 
the projections. We are interested in the problem of reconstructing a tiling which 
has given projections. Some simple variants of this problem, involving tiles that are 
(f-) . 1 x 1 or 1 x 2 rectangles, have been studied in the past, and were proved to be either 

solvable in polynomial time or NP-complete. In this note we make progress toward 
a comprehensive classification of various tiling reconstruction problems, by proving 
. NP-completeness results for several sets of tiles. 

00 



1 Introduction 



CZ3 

In Discrete Tomography we want to reconstruct a discrete object from its projections. This 
paper is concerned with the reconstruction of tilings. We are given a collection of tiles, 
where each tile can have a different shape. A tiling is a placement of non-overlapping copies 
of the tiles on a n x n grid, where each copy is obtained by translating one of the tiles. (In 
this note we do not allow tile rotations, although one could also consider the variant with 
rotations.) Projections of a tiling determine the number of tiles intersected by each row 
and column. Given such projections, we wish to reconstruct a tiling consistent with these 
projections, or to report that such a tiling does not exist. 

Formally, a tile t is denned to be a finite subset of Z 2 . In this paper we only consider 
tiles that are hole-less polyominoes. By + 1 = {(i + i' ,j + j') : (i',f) G t] we denote 
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the translation of t by vector G Z 2 . Fix a finite multiset of tiles T = {ti, £2, ■ ■ ■ , th}- 
Without loss of generality we assume that every tile tk contains (0, 0), the so-called center 
of the tile, and in this paper it will always be the upper-left corner. We refer to the index 
k as the type of the tile. The tiles are identified by their type, and different tiles may have 
the same shape. One can think of tiles which are of the same shape but of different types 
as being of different colors. 
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Figure 1: A tiling of the 8x8 grid with its projections. The tile centers are marked by a 
circle. 

A T -tiling of a grid G = Z n x Z n is a finite set T C Z n x Z„ x [l,h], such that 
the sets + t k , for all (i,j,k) G T, are disjoint and contained in G. If T = {t}, we 
will sometimes write simply t-tiling instead of {£}-tiling. The center projections of T are 
vectors r, c G M nx , where 

n,k = k) G T}\ and 

c j,k = \{(i,j,k) G T}\. 

The numbers count the number of tiles of type k whose center is in row i, the numbers 
cj 7 k count the same for column j. In a similar manner we define the cell projections r', d of 
T, where we count for each row, each column, and each type the number of cells covered 
by that type of tile. 

If tilings T and T' are disjoint, then projections of T U T' are the sums of the pro- 
jections of T and of T' . (This is true for both types of projections.) Therefore the set 
of projections of all tilings T with a single tile (|T| = 1) spans the set of all projections. 
The canonical bijection between single-tile center projections and single-tile cell projections 
implies a bijection between all center and cell projections. From now on we will use the 
term "projection" for center projections, unless stated explicitly otherwise. 

Note that we do not require the tilings to cover the whole grid. Tilings that cover the 
whole grid are called complete. Each tiling problem can be mapped into an equivalent 
complete tiling by adding one "clear" lxl square tile whose row cell projections are n 
minus the total sums of the other tiles' row cell projections, and the column projections 
are defined analogously. 

Figure [I] illustrates this definition for the tile set T = {£i, ^2}, where £1 is the 3x3 
square and t<i is the 2x2 square. The numbers on the left are the projections and the 
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numbers on the top are the Cj t k- Columns are numbered from left to right and rows from 
top to bottom, with indices ranging from to 7. For example c^i = 2 because column 4 
contains two centers of tile t\. 

Given a tiling T, the computation of its projections is straightforward. Consider now 
the inverse problem: given the vectors r, c, find a T-tiling T with projections r, c. This is 
called the T -reconstruction problem. The related decision problem ("is there such a tiling 
T?") is called the T -consistency problem, or simply the T-tiling problem. 



Our results. For some types of tiles the reconstruction problem is easy to solve, while for 
other it may be hard. Table |l] summarizes the complexity of various tiling reconstruction 
problems, including both our results and previous work. In this table, by "NP-complete" 
we mean that the consistency problem is NP-complete. 



The /-atom reconstruction problem. The simplest tile is a 1 x 1 square, which we call 
a cell or an atom (the original motivation for this problem came from the reconstruction 
of polyatomic structures). When T consists of I different cells, we will refer to the T-tiling 
problem as the l-atom problem. Reconstructing 1-atom tilings is easy and can be solved 
in time linear in the size of the grid, as shown in 1957 by Ryser [|Rys57| , [Rys63|| . For 3 
or more atoms (cells of different type), the reconstruction problem is NP-hard | CD01|| (see 
also [|CCP00||). For 2 atoms, the complexity of the problem remains open. 



One tile. For a single tile, it is known that if the tile is a horizontal bar, i.e. a rectangle of 
height one, the problem is as easy as reconstructing 1-atom tilings [ UCRK , PicOl j. There 
exist other types of tiles, however, for which the problem is NP-complete. Two such tiles 
are given in Table [l]. The problem remains open for rectangular tiles, even for the 2x2 
square. 



Two tiles. For pairs of tiles the situation is quite different. For horizontal and vertical 
dominoes — 1x2 and 2x1 rectangles — the problem is NP-hard. However, the case when 
the domino tiling is required to be complete is open. The problem is also open for vertical 
dominoes and single cells. For squares the problem is NP-hard, both for two types of 2 x 2 
squares, and for a single 2x2 square and a cell. 



2 NP-hardness Proofs 

We now present our NP-hardness results. In our proofs we reduce the 3-atom problem to 
the given version of the T-tiling problem. A similar strategy was used earlier by Diirr, 
Goles, Rapaport and Remila |pGRR|| to show that reconstructing tilings of given sub-grids 



- grids with forbidden regions — with only vertical and horizontal dominoes is NP-hard 
(even if the tilings are required to be complete). 

The general idea of the proofs can be summarized as follows. We think about the 
3-atom problem as a 4-atom complete tiling problem, with an additional "clear" atom. 
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Table 1: The complexity of different versions of the reconstruction/consistency problem. 



For convenience, we name each possible tile in the 3-atom problem as yellow, blue, red 
or clear. Throughout this section, by (r, c) we will denote the given instance of the 3- 
atom problem. We will map (r, c) into an instance (r', d) of the T-tiling problem under 
consideration. In all proofs we assume, without loss of generality, that r ijfc = V. c^k 
for all k. This assumption is valid, since we can extend any mapping to instances in which 
Yli r i,k 7^ Ylj c j,k> by mapping them into an arbitrary fixed negative instance (r',d), say 
to one in which the totals of row sums are not equal to the totals of column sums. This 
does not affect the asymptotic running time nor the correctness of the transformation. 

For simplicity, assume first that T contains just one tile. To construct (r',d), we 
choose a small dxd grid B, called a block, that can be tiled in only four possible ways (this 
restriction will be relaxed in some proofs). Each of these four so-called admissible block 
tilings will correspond to one atom. Instance (r', d) will have grid dimensions nd x nd. 
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We view this grid as an n x n matrix consisting of d x d blocks. A segment of rows 
numbered id, . . . , (i + l)d — 1 will be referred to as a block-row i. The transformation 
maps r^i, . . . , 7^4, that is, the atom projections of row i, into a length-a? vector which is a 
projection of block-row i. This vector is the linear combination of the projections of the 
admissible tilings of B with coefficients r^i, . . . , r^. The column projections are mapped 
in the same way. 

If T has h > 1 tiles, the transformation is the same, except that now the block projec- 
tions are not length-d vectors but d x h matrices. 

Obviously, for any tile set T, the T-consistency problem is in NP. For any choice of B 
and its admissible tilings, the method outlined above can be implemented in polynomial 
time. It also has the property that if (r, c) has a solution then so does {r',d). For if T 
is a solution of the 3-atom problem with projections (r, c), then the tiling T" obtained by 
replacing each atom by its corresponding admissible block is a T-tiling with projections 
(r',d). Thus the above ingredients of NP-completeness arguments will be omitted in the 
proofs below, and we will focus exclusively on proving the following implication: if (r', d) 
has a solution then (r, c) has a solution. 

The main difficulty is to construct B to make this latest implication work. In other 
words, we need the property that any tiling of the resulting instance (r', d) can be trans- 
formed back into a solution of (r,c). To achieve this, we choose B and the admissible 
tilings so that the following two conditions hold: 

(npcl) any tiling with projections (r',d) consists only of admissible block tilings, and 

(npc2) from the projections of the block-rows we can uniquely extract the projections of 
the atoms in the corresponding rows of the 3-atom problem. 

To enforce condition (npcl), we use techniques inspired by classical structure theorems 
for realizations of — 1 matrices with given projections [ Rys63[ . Another useful method 



involves the total counts of different colors. By Y, B, R and C we will denote the total 
number of yellow, blue, red and clear atoms in (r, c). We have B + Y + R + C = n 2 . The 
block projections impose additional restrictions on how many centers of the tiles in T can 
occur on certain positions in the blocks. These restrictions can be expressed in terms of 
numbers Y, B, R and C. By investigating these constraints, we prove that non-admissible 
tilings cannot occur. 

We now discuss condition (npc2). Suppose that there is a T-tiling X" with projections 
(r', d). By (npcl), each block in T' is admissible. We transform X" into a solution T of the 
3-atom problem by replacing each admissible block by its corresponding atom. To satisfy 
(npc2), we need to show that the projections of T are (r, c). 

Number the admissible tilings from 1 to 4 and name their row projections e l5 e 2 , e 3 , e 4 . 
Let hi be the projection of block-row i in T' and qj the number of blocks in block-row i 
with the jth admissible tiling. Then the numbers qj satisfy: 

giei + q 2 e 2 + <? 3 e 3 + <?4e 4 = bj. (1) 



5 



By the construction, equation (|1|) has a solution qj = r it j, for j = 1,2,3,4. For (npc2) 
to hold, we need to ensure that ([!]) does not have any other solutions in non-negative 
integers. This can be easily accomplished by choosing the admissible tilings for which the 
projections ei,e 2 ,e 3 ,e4 are linearly independent. 

In fact (npc2) will hold even for a weaker condition. Note that the numbers satisfy 
r i,i + r i,2 + r i,3 + r «,4 = n - So we extend each e 3 - by adding to it one coordinate with value 
1, and we similarly extend each vector bj by adding to it one coordinate with value n. 
(Technically, the and bj can be d x h matrices, but for the purpose of the transformation 
we can as well treat them as vectors of length dh. Then the extended vectors will have 
length dh + 1.) If these new vectors ei, e2, e3, are linearly independent, we can use our 
admissible tilings for the transformation. Although we do not use it in the paper, it may 
be worth to mention that the linear independence of these extended vectors is equivalent 
to a condition called affine linear independence (see page 3 of ||GLS88|| ). 



[o 



Theorem 1 Let t be the tile 




The consistency problem for t is MP- complete. 



Proof. The proof follows the method outlined above. We reduce the 3-atom consistency 
problem to the t-tiling. We treat the 3-atom problems as the (equivalent) complete 4-atom 
problem, by adding an extra "clear" cell tile. We use a block of size 7x7. The admissible 
tilings are the tilings of the 7x7 block with three tiles t. There are exactly four admissible 
tilings. These tilings and their associations to different atoms are shown in Figure |2]. Using 
the projections of these admissible block tilings, we map any instance r, c G N nx4 of the 
3-atom consistency problem into an instance r',c' G N 7n of the t-tiling problem. 
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Figure 2: The four admissible tilings of a block with 3 tiles t. 

The transformation works as follows. Let ei,e2,e3,e4 G N 7 denote the row projection 
vectors of the four tilings in Figure 0. For every row i, its projections form a 4-dimensional 
vector (^1,^2,^3,^4). We will map it into N 7 . The projections r' in the resulting 
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instance are defined by 



7i 



7i+l 
i 

7i+2 
I 

7i+3 

r 7i+4 
l 
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7i+5 



J"»,i e i + r i^2 + ri 3 e 3 + r iA e 4 . 



\ r 7i +e J 



In other words we set r' 7i 



r i,l + n i r 7i+l — r i,2i r 7i+2 — r i,3; 



7i+3 



iA 



+ n and r'-. 



7i+l 







for I = 4, 5, 6. The column projections d are determined in a similar manner. 

We need to show that if (r', d) has a solution then (r, c) also has a solution. To this 
end, let T be an arbitrary solution to (r',d). We claim that in T every block is in one 
of the four configurations of Figure 0. This is true because the rows and columns whose 
indices modulo 7 are greater than 3 have projection and therefore all tiles t are strictly 
contained in a block. Further, (r f , d) requires 3n 2 tiles t in total, and each of the n 2 blocks 
contains at most three tiles of type t. 
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Figure 3: Reducing the 3-atom reconstruction problem to the t-tiling reconstruction prob- 
lem. 

The vectors are linearly independent. So the projections (r' 7i+k )kez 7 of a block row i 
can be uniquely written as Ylt=i r «,fc e fe- This ensures that T can be mapped into a 3-atom 
tiling with projections r, c (see Figure [3]). □ 

For the next result we need the following classical result on the structure of — 1 



matrices with given projections [ |Rys63|| . We state this fact in terms of tiling with cells. 



Fact 1 Let r, c 6 Z n be an n x n instance of the tiling problem with cells. Let I C Z n be 
a row set and J C Z n a column set. If 



J2n-J2 c 3 = \ I x J \, 



then in every solution the set I x J must be completely tiled and the set I x J must be 
completely empty. 
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Proof. Let T be a solution to (r, c). Let a be the number of cells in / x J, b the number 
of cells in J x J and c the number of cells in / x J. Then 

rj — Cj = (a + fo) — (6 + c) = a — c. 
If a — c = \I x J| then a = |J x J\ and c = 0, which concludes the proof. □ 
Theorem 2 The consistency problem for 1x2 and 2x1 rectangles (dominoes) is NP- 

Proof. We follow the idea outlined at the beginning of this section. We use a 3 x 3 block. 
The four admissible tilings of the block are shown in Figure |j. 
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Figure 4: Four admissible tilings of the 3x3 grid with dominoes. 

We need to show that if (r ; , d) has a solution then (r, c) has a solution. Note that the row 
and column projection matrices of the four tilings in Figure (|are not linearly independent, 
but at least satisfy the weaker condition described on page |5], which is enough for the 
reduction. Therefore to complete the proof we need to show that in any solution of the 
tiling instance every block is admissible. 

Let I = J = {i G Z 3n : i mod 3 > 0} be row and column sets. Denote the "clear" cell 
by t 3 . Recall that Y,B,R and C denote the total numbers of yellow, blue, red and clear 
atoms in (r,c). Then we have r' i 3 - £\ £ 7 c' j 3 = (5Y + 5B + 4R+6C) -(Y + B + 2C) = 
AY + AB + AR + AC = An 2 . Fact | implies that in every block dominoes can only appear 
in the first row or first column of each block and the top-left cell is always covered by a 
domino. All the tilings that satisfy these condition are the admissible tilings in Figure |]. 

□ 

Theorem 3 The consistency problem for two types of 2 x 2 squares is NF-complete. 

Proof. We refer to the two types of 2 x 2 squares as light and dark squares. We use the 
4x4 block and four admissible block tilings shown in Figure |5|. The row and cell projection 
matrices of the admissible block tilings are linearly independent. Thus, to complete the 
proof, we have to show that any solution of (r', d) uses only the four admissible block 
tilings. 
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Figure 5: Four tilings of the 4x4 grid with two types of squares. 



If we consider the cell projections rather than the center projections, we see that row 
1 and column 1 of every block must be completely tiled with light squares (recall that 
rows and columns are numbered from 0). Therefore a block can only be in one of the five 
tilings shown in Figure j|. The fifth tiling — which is called bad tiling — has the same 
row projections as a "yellow" block and the same column projections as a "blue" block. 
By column projections for columns j = 1 (mod 4), in any tiling there will be Y "yellow" 
blocks, and by row projections for rows i — 1 (mod 4), there will be B "blue" blocks. There 
are C + R remaining blocks, and 4(C + R) tiles must appear in these remaining blocks, so 
each of these remaining blocks must have four tiles. Thus the bad tiling cannot occur. □ 

The same technique can be used to reduce the 2-atom problem to the single-type 
square tiling problem. In this reduction, only the first 3 block tilings of Figure |5] are used. 
It can also be used to prove MP-completeness of the cell-and-square tiling problem. In 
this reduction, the dark square is replaced by a cell in Figure [|, without modifying the 
projections. 

Theorem 4 If the consistency problem for 2x2 squares can be solved in polynomial time, 
then the 2-atom problem can be solved in polynomial time. 



Theorem 5 The consistency problem for lxl cells and 2x2 squares is NF-complete. 



Theorem 6 Let t be the L-shaped tile 



The consistency problem for t is NF-complete. 



Proof. We reduce the 3-atom problem to the t-tiling problem using the admissible block 
tilings of Figure ^. The first three tilings correspond to the colored atoms, and the two 
last tilings (with identical projections) correspond to the clear atom. The row and column 
projections of the admissible tilings are linearly independent. 

It is sufficient to show that in any solution T to (r 1 , d) every block must be admissible. 
We define the matrix M G N 6x6 where m it j is the number of tiles in T whose center equals 
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Figure 6: Five admissible tilings of the 6x6 block with the L-shaped tile. The last two 
tilings have the same projections. 



modulo 6. Row and column sums of T imply that M must have the following form 
(the numbers on the left and on top are the row and column sums of M): 
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The tiles centered at and (i — 1, j — 1) overlap, and they both overlap with the tile 
centered at (i — 1, j) or (i,j — 1). This introduces two additional constraints on the matrix: 
for every i, j 

m i-i,j-i + m i-ij + m i,j — n 2 ? an d (2) 
mi_i,j_i + rriij-t + m itj < n 2 , (3) 

where we set rriij = for i = — 1 or j = — 1. 

Every block must have 2 centers in row 4 (recall that rows and columns are numbered 
from 0). So row 3 can have at most 1 center. We now consider sub-blocks that consist of 
rows 3,4, 5. By the above argument, the tiles that are fully contained in these sub-blocks 
must have one of the following configurations: 





2 




3 




4 






5 














o]_ 




o]_ 


o 




o 































This immediately gives 7713,0 = 7773,2 = 0, 7774,0 — n. Let a s be the number of blocks 
whose last three rows are in configuration of type s above, for s — 1, 5. Then 0,4, = 7773,3 
and a 2 + a 3 + a 5 = 7774 3. Since the projections of row 3 and column 3 of M are equal C, 
we get a 4 + a 5 = C and a 2 + 03 + a 4 + a 5 + 777 0i3 + 777.1,3 + m 2,3 = C. Thus a 2 = a 3 = m 0j3 = 
777i,3 — m 2,3 — 0. This means that sub-blocks of types 2 and 3 cannot occur and that M 
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has the form 
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where we write a = as, for simplicity. 

Projections of row and column imply mo,o < R and m^o + "ii,o + m 2,o = n 2 + R, 
and from fl2|) we have m^o + m.2,Q < ?i 2 - Therefore mo,o = -R and m^o + 777.2,0 = ti 2 - The 
first equation forces m ,i = m , 2 = 0, while the second forces m^i = 0, by fl3|). Hence 
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for some integer b > 0. Projections of column 1 and rows 1, 2 imply that 

( R 
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By inequality (Q) for (hj) — (3, 1) we have n 2 — b + a < -a 2 , so a = 6, because all entries 
are non-negative. Thus m 12 = 0. 

Write >c (i',f) if each block has a center at exactly one of the positions or 
(i',f)- Clearly, if rriij + rrii'j/ = n 2 and the tiles centered at (i',f) overlap, then 

cxi Since a = 6, by the form of M above, we get the following relations: 

(1,0) tx (2,0) M (3,1) X (4,2) txl (4,3) and (1,1) M (2,2). 

Write = if each block has a tile centered at iff it has a tile centered at 



(i By the above, we get 

(1,0) = (3,1) = (4,3) and (2,0) 



(4,2) 



By extending the three allowed configurations (number 1,4,5) of the rows 3, 4, 5, and using 
the above constraints, we obtain that the only block tilings that meet these conditions are 
the admissible tilings in Figure E. □ 
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3 Conclusion 



We proved NP-completeness for several variants of the tiling problem, but what can be 
said about the infinitely many variants for which the complexity remains open? 

A tile t of width w and height h is said to be interlocking if there is a box of width < 2w 
and of height < 2h which contains two disjoint copies of t. For example, the tiles from 
Theorems [I] and |6] are interlocking, while rectangles are not. There are non-rectangular 
tiles that are not interlocking, for example the U-shaped tiles. For instances consisting of 
one interlocking tile we believe the tiling problem to be NP-complete. 

We are also confident that the problem is NP-complete for all variants involving at least 
two different tiles, one of width > 2 and one (possibly the same one) of height > 2. This 
condition ensures that the problem is not invariant under column or row permutations. 

We believe that the techniques developed in this paper will be useful in developing 
generic transformations that can show NP-completeness of wide classes of tiling problems. 
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